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Abstract 



We report a determination of the B^-B^ mixing parameter Amj based on 
the time evolution of dilepton yields in T(45) decays. The measurement is 
based on a 5.9 fb _1 data sample collected by the Belle detector at KEKB. The 
proper-time difference distributions for same-sign and opposite-sign dilepton 
events are simultaneously fitted to an expression containing Am^ as a free 
parameter. Using both muons and electrons, we obtain Am^ = 0.463 ± 
0.008 (stat.) ±0.016 (sys.) ps~ 4 . This is the first determination of Amd from 
time evolution measurements at the T(45). We also place limits on possible 
CPT violations. 

PACS numbers: 13.20.H 
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Typeset using REVTeX 



The frequency of B®-B® mixing is proportional to the mass difference between the two 
mass eigenstates of the neutral B meson, Amd, and is a fundamental parameter of the B 
system. Measurements of Am^ derived from the time evolution of B® decays have been 
reported by CDF, SLD, and the LEP experiments ][]; ARGUS and CLEO have measured it 
using the integrated fraction of same-flavor B pair decays in T(4S) events We report 

here the first determination of Am^ based on the time evolution of B® decays in T(4S f ) 
events produced in asymmetric e + e~ collisions, using data collected by the Belle detector 
H at the KEKB storage ring [|. 

At the T(45), the asymmetry in time evolution between same-flavor (B®B®, B^B®) and 
opposite-flavor (B®B®) decay pairs exhibits an oscillation as a function of the proper time 
difference between the two 5-meson decays, At, with a frequency that is proportional to 
Amj. In KEKB, collisions between 8.0 GeV electrons and 3.5 GeV positrons have a center 
of mass (CM) motion along the electron beam direction (z direction) with a Lorentz boost of 
7/3 = 0.425. Since each of the two B's is produced nearly at rest in the CM, the separation 
of their decay vertices in the lab frame is proportional to At and has an average magnitude 
of 200/iin. High-momentum leptons can be used both for tagging the B flavor and for 
determining the decay vertex with good accuracy. The At in dilepton events can thus be 
used to measure the time evolution of B decays. The same analysis can be used to test CPT 
conservation by the inclusion of the complex parameter cos# in the fit 0. 

The analysis presented here is based on integrated luminosities of 5.9 fb^ 1 at the T(45) 
resonance and 0.6 fb _1 at an energy that is 60 MeV below the peak. 

The Belle detector consists of a silicon vertex detector (SVD), a central drift cham- 
ber (CDC), an array of 1188 aerogel Cerenkov counters (ACC), 128 time-of-flight (TOF) 
scintillation counters, and an electromagnetic calorimeter containing 8736 CsI(Tl) crystals 
(ECL), all located inside the 3.4-m-diameter superconducting solenoid that generates a 1.5 T 
magnetic field. An iron return yoke, outside the solenoid, is segmented into 14 layers of 4.7- 
cm-thick iron plates alternating with a system of resistive plate counters that is used to 
identify muons and Ki mesons (KLM). 

Hadronic events are required to have at least five tracks, an event vertex with radial 
and z coordinates respectively within 1.5 cm and 3.5 cm of the origin, a total reconstructed 
CM energy greater than 0.5W 7 (W is the T(AS) CM energy), a ^-component of the net 
reconstructed CM momentum less than 0.3W/c, a total CM calorimeter energy between 
0.0251^ and 0.9011 7 , and a ratio Ri of the second and zeroth Fox- Wolfram moments |7j that 
is less than 0.7. While the R2 cut suppresses events of non-T(4S f ) origin, all other cuts are 
intended to remove the beam-related background and QED events. 

For electron identification, we use position, cluster energy, and shower shape in the ECL, 
dE/dx in the CDC, and hit information in the ACC. This is ~ 90% efficient for electrons 
and has a ~ 0.3% misidentification probability for charged hadrons with momenta above 1 
GeV/c. Electrons from 7 conversions are removed. 

Muon selection is based on KLM hits associated with charged tracks. The range of the 
tracks and the matching quality of the hits are used. The efficiency is ~ 85% for muons 
with momentum above 1 GeV/c and the misidentification probability is ~ 2%. 

Events containing leptons from J/ip decays are rejected. In addition, lepton candidates 
are required to satisfy 30° < 6 < 135°; 1.1 GeV/c < p* < 2.3 GeV/c; | dr 1F |< 0.05 cm; 
I dz lp |< 2.0 cm; and have at least one (two) associated SVD hit(s) in the r-<p (r-z) view, 
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where 6 is the laboratory polar angle, p* is the CM momentum, and dr IP and dz lp are the 
distances of closest approach to the run-dependent interaction point. To reduce secondary 
leptons and fakes from the same B and from the continuum, which tend to be back-to- 
back, the opening angle 9* u between the leptons in the CM frame is required to satisfy 
—0.8 < cos#££ < 0.95. The application of the above-listed criteria yields 8573 same-sign 
(SS) and 40981 opposite-sign (OS) dilepton events on the T(4S f ), and 40 SS and 198 OS 
dilepton events below the resonance. 

The z-vertex of leptons is determined from the intersection of the lepton tracks with 
the profile of B® decay vertices, which is estimated from the profile of the beam interaction 
point (IP) convolved with the average B flight length (~ 20/im in the T(4S I ) rest frame). 
The mean position and the width (cr^, 0"^ p , of the IP are determined on a run-by-run 
basis using hadronic events. We find cx^ p = 100-120/mi, <r^ p ~ 5/im and a 1 ^ = 2-3mm. The 
proper-time difference is calculated from the z positions of the two lepton vertices using the 
relation At = /S.z/cf3 r j ) where Az = z\ — z<i is the difference between the two z vertices. 
For OS events, the positively charged lepton is taken as the first lepton {z\). For SS events 
the absolute value of Az is used. 

The observed SS and OS dilepton proper-time distributions have contributions from 
"signal," defined as events where both leptons are primary leptons from semileptonic decay 
of B\ or B + , and "background," where at least one lepton is secondary or fake, or the 
event is from the non-T(4S') continuum. The value of Arrid was extracted by simultaneously 
fitting the two distributions to the respective sums of contributions from all known signal 
and background sources. 

Each dilepton is identified with one of the event types listed in Table |. Each event type 
is categorized as either signal (S), correctly tagged background (C), or incorrectly tagged 
background (W). For each, we parameterize the proper-time distribution as the product of 
the number of contributing events (N) , an overall selection efficiency (e) , and a normalized 
distribution function P(At). The values of N depend on the total number of T(4S) events in 
the data sample (N^s), the branching fractions of the T(4S) to neutral and charged B pairs 
(fo and /+ = 1 — /o) , the semileptonic branching fractions (&o and b + for neutral and charged 
B), and, for neutral S's, the mixed event fraction \d = x 2 d / (2(1 + x 2 d )) where Xd = T B oArrid 
and r B o is the B d lifetime. The evaluations of iV are summarized in Table |. The efficiencies 
e are determined by Monte Carlo (MC) simulation. 

The observed proper-time distribution function P(At) for the signal is a convolution of 
a root distribution with a detector response function, g, and is given by: 

p(Af\ = f d(At - At')F(At')d(At') 
1 ' f J g(At — At')F(At')d(At')d(At) ' { ) 

where the respective theoretical root functions for mixed B°, unmixed B°, and charged B 
are 

F(At) = (l/4r B o)e HA<l/rs d [1 - cos(Am d At)} (2) 

F(At) = (l/4r B o)e~' A ' l/Ts d [1 + cos(Am d At)] (3) 

F(At) = (l/2r B+ )e- |A * l/T a+. (4) 

MC simulations of generic B + B~, unmixed and mixed B°, and continuum events are 
used to determine the background Az distributions. The dominant background source is a 
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primary lepton paired with a secondary lepton from a c-quark. The shape as well as the 
normalization of the background from neutral B events depends on Am^. To account for 
this, we generated two samples of generic neutral B events, one with Arrid = 0.464 ps _1 and 
one with Am<j = 0.423 ps _1 . Background distributions for arbitrary Am^ are determined 
by linear interpolation. 

We use the Az distribution of dileptons from J /if) decays in the data for g; for these 
events the root distribution is a delta function and the lepton momentum spectra are in the 
same region as those of primary leptons from B decays. The Az distribution for J /if) events, 
which has a = 112 /xm, agrees with the MC distribution if it is convolved with a Gaussian of 
a = 50 ± 18 /im. This is due to an imperfect detector simulation, and we correct this effect 
by applying a convolution with o = 50 /im to each MC-determined background distribution. 

To extract Am^, a binned maximum likelihood fit is performed simultaneously to the 
Az distributions of the SS and OS dileptons. Each fitting function is a sum of signal 
and background distributions. In order to properly take into account the tails of the Az 
distributions, the signal response function and the background distribution are given in the 
form of a lookup table rather than an analytic function. We fix the parameters t b o = 1.548 
ps 0, /+//„ = 1.05 0, and t b +/t b o = 1.06 0, and limit the fit region to |A*| < 1.85 mm. 
The constraint b + /bo = t b +/t b o is imposed. The continuum contribution is fixed to that 
of off-resonance data, scaled to account for luminosity and energy differences. The relative 
selection efficiencies for the event types (mixed B%, unmixed B®, and charged B) within 
each tag type (S, C, and W) are fixed, resulting in two free parameters (efficiency ratios 
in C/S and W/S) in addition to Arrid and the overall normalization. The fit result is 
Am d = 0.463 ± 0.008 ps" 1 with x 2 /DOF = 333/376. The efficiency ratios in C/S and W/S 
are (9.66 ± 1.39) x 10~ 3 and (6.98 ± 0.25) x 10~ 3 , respectively, which give signal fractions 
to be 32.1%(SS) and 77.5%(OS). Figure [l] shows the Az distributions for the data together 
with the fitted curves. Fig. ^| shows the OS and SS asymmetry, (N Q s — N ss )/(N s + N ss ), 
for data together with the result of the fit. 

As a cross check, we also measured Am^ using a fitting method that differs from the one 



described above in the following aspects [ 10] : a) an unbinned rather than binned maximum 
likelihood fit; b) response function is the sum of three Gaussians, with parameters determined 
from the dileptons from J /if) decays; c) backgrounds separated into SS, OS, rather than C, 
W; d) background distributions were analytic functions, with parameters determined by 
fitting to MC We find Am^ in the range 0.460 to 0.483 ps -1 depending on the choice of 
analytic forms for the backgrounds, which is consistent with the primary result. 

The systematic errors were estimated by repeating the fits for different input parameters. 
The main contribution originates from uncertainties on input parameters and from deter- 
mination of the response functions. Contributions of /+/ f Q , r B o and r B +/r B o are estimated 

d d 

by adjusting each in turn by the amount of its uncertainty. Contribution of the response 
function arises from the possibility that it differs from the true dilepton response function, 
from the statistical uncertainties of the determination, and from the fact that the calculation 
of proper time At = Az/cf3^ is not exact due to the motion of the B's in the CM and 
the energy spread of the beams. To estimate the first possibility, we used the MC dilepton 
response function convolved with a Gaussian of a = 50/im. For the second, we varied the 
number of entries on a bin-by-bin basis by the amount of the statistical errors. For the third, 
we compared two fits, one using a response function obtained for the true At difference and 
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a second obtained for Az. 

We also consider the uncertainty from the background simulation. We assigned a 35% 
error for the fake rate and adjusted the fake rate by this amount. We varied the branching 
ratios of B decaying to D° and D + in the MC in accord with the experimental uncertainties 
|Tl"[ . We varied the width of the Gaussian used to correct for an imperfect detector simulation 



by ±18/im. 

It is assumed in the fit that Ar, the difference between the decay widths of the neutral 
B mass eigenstates, is zero. Although no significant experimental constraint exists ||, it is 
predicted based on solid theoretical grounds to be very small {AT/T < 1%) ||. We repeated 
the fit including the effects of AT/T = 1% and found the shift in the result to be negligible 
(<0.001 ps- 1 ). 

Contributions to the systematic error from the above sources are summarized in Table O. 
The total systematic error is obtained by summing all errors in quadrature: 

Am d = 0.463 ± 0.008 (stat) ± 0.016 (sys) ps -1 . 

When the constraint of CPT conservation is removed in B^-B^ mixing, the theoretical 
functions (§) and (§) are modified and become 

F(At) = (| sin#|74r B o)e~ |A<l/Ts d[l _ C0S {Am d At)] (5) 
F(At) = (l/4T B o)e HA * l/r ^[l + | cos#| 2 



d ' 

|2\ 



+(1 - I cosfln cos(Am d At) 
-2Im(cos0)sin(Am d At)] (6) 

and Xd becomes (| sin#| 2 a^)/[| sin9\ 2 x d + (2 + x\ + x 2 d \ cos^l 2 )]. A non-zero value of the 
complex parameter cos# would be an indication of CPT violation. The result of the fit is 

m 



Im{cos9) = 0.035 ± 0.029(stat) ± 0.051(sys) 
Re{cos6) = 0.00 ± 0.15(stat) ± 0.06(sys) 

and Amd = 0.461 ps -1 . These results imply || the upper limits \m B o — m B o\/m B o < 
1.6 x 10~ 14 and \Y B0 - T B o\/T B o < 0.161 at the 90% C.L. 

In summary, we report the first determination of Airtd using the time evolution of B°- 
mesons produced in T(45) decays. We obtain Am<j = 0.463±0.008 (stat) ±0.016 (sys) ps -1 , 
which is consistent with the world average value Amd — 0.472 ± 0.017 ps -1 [§. We have 
also examined CPT violation and obtain the first limit on (m B o — m B o)/m B o and a limit on 
(r#o — r^ )/rso that is compatible with the previous measurement [|13[ . 
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TABLE I. Categorization of event types contributiong to dilepton events. 



11 event type 




A 


tag type 


SS signal 


B^, mixed 




s 


SS background 


B%, mixed 


NisfoXd 


c 




7-lfl - 1 

Bj, unmixed 


N AS fo{l ~ Xd) 


w 




B + B~ 


N 4S f+ 


w 




continuum 






OS signal 


B®, unmixed 


JV4s/o(l " Xd)bl 


s 




B+B- 




s 


OS background 


B^, mixed 


N 4 sfoXd 


w 




B^, unmixed 


N 4S fo(l ~ Xd) 


c 




B+B- 


N 4 sU 


c 




continuum 


A con t 





TABLE II. Summary of systematic errors. 



Source (uncertainty) 


Am d 


/+//„ (± 0.08) 


±0.009 


B° d lifetime (± 0.032 ps) 


±0.004 


t b+ /t b o (± 0.03) 

a 


±0.009 


response function 


±0.005 


background fake rate (±35%) 


±0.004 


B{B -» DX) (D°: ±4.6%, D+: ±14.3%) 


±0.002 


continuum (SS: ±16%, OS: ±7%) 


±0.002 


detector resolution, BG (±18 jum) 


±0.001 


Monte Carlo statistics 


±0.004 


total 


±0.016 
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FIG. 1. Az distribution of dileptons for data together with fit result. The upper plot shows the 
distributions for same-sign, and the lower plot for opposite-sign dileptons. Signal and background 
dileptons obtained from the fit are also shown. 
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FIG. 2. Opposite and same-sign dilepton asymmetry vs Az. The asymmetry is defined as 
A(Az) = (Nqs — ^* r ss)/(-^OS + N$s). The points are the data. The curve is the result of the fit. 
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